(3.K  IIBM* 


awnsHWis 


The  University  of  Alberto 
Printing  Department 
Edmonton,  Alberta 


■ 


THE  UNIVERSITY  OF  ALBERTA 

LIAPUNOV  STABILITY  OF  DISTRIBUTED 

PARAMETER  SYSTEMS 


by 


© 


Thomas  Wayne 


Lee 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  MASTER  OF  SCIENCE 


DEPARTMENT  OF  ELECTRICAL  ENGINEERING 


EDMONTON,  ALBERTA 
FALL,  1970 


UNIVERSITY  OF  ALBERTA 
FACULTY  OF  GRADUATE  STUDIES 


The  undersigned  certify  that  they  have 
read,  and  recommend  to  the  Faculty  of  Graduate 
Studies  for  acceptance,  a  thesis  entitled  Liapunov 
Stability  of  Distributed  Parameter  Systems  sub¬ 
mitted  by  T.  Wayne  Lee  in  partial  fulfilment  of 
the  requirements  for  the  degree  of  Master  of 
Science. 


T\ 


ABSTRACT 


During  the  past  few  years,  interest  in  studying 
stability  of  distributed  parameter  systems  via  Liapunov's 
second  method  has  slowly  been  forming.  In  this  thesis  two 
techniques  are  developed  for  finding  Liapunov  functionals 
for  certain  types  of  distributed  systems.  In  the  first 
technique,  the  use  of  a  non-identity  matrix  in  a  Liapunov 
functional  of  quadratic  form  proves  to  be  more  successful 
than  methods  used  in  which  the  weighting  matrix  was 
replaced  by  the  identity  matrix.  The  second  technique  is 
an  extension  to  the  variable  gradient  method  of  generating 
Liapunov  functions.  A  systematic  algorithm  analogous  to 
the  one  used  for  lumped  systems  is  developed  and  illustrated 
by  two  examples. 
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CHAPTER  1 

REVIEW  OF  STABILITY  METHODS 
FOR  DISTRIBUTED  PARAMETER  SYSTEMS 

1.1  Need  For  Stability  Analysis 

A  system  is  said  to  be  a  distributed  parameter 
system  if  its  variables  depend  both  on  time  and  the  spatial 
coordinates.  If  the  spatial  dependence  is  sufficiently 
small  that  it  can  be  neglected,  then  the  system  is  called 
a  lumped  parameter  system.  For  a  distributed  parameter 
system,  usually  the  governing  equations  are  partial 
differential  equations,  while,  for  a  lumped  system,  these 
equations  reduce  to  ordinary  differential  equations. 

One  of  the  difficulties  in  controlling  physical 
systems  is  encountered  in  assuring  adequate  stability.  A 
qualitative  definition  of  stability  is  as  follows 
DEFINITION 

A  system  with  a  suitable  reponse  to  a  set  of 

inputs  and  initial  conditions,  is  said  to  be  stable  if  it 

returns  to  a  new  response,  close  to  the  original,  whenever 

small  changes  occur  in  the  inputs  and/or  initial  conditions. 

Currently,  much  of  the  work  on  stability  of 

distributed  parameter  systems  comes  from  the  area  of  process 

1  *  •  19 

control.  The  review  papers  of  Aris  and  Lapidus  cover 


Numbers  placed  above  the  line  of  text  refer  to  the  references. 
Only  those  works  relevant  to  this  thesis  are  listed.  For  an 
extensive  bibliography  on  stability  of  distributed  parameter 
systems  see  Wang^3 . 
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the  major  works,  many  of  which  were  published  within  the 

past  five  years.  Systems  in  process  control  which  are 

characterized  by,  say,  diffusion  or  conduction  are  well 

suited  to  be  analysed  as  distributed  parameter  systems. 

The  need  for  stability  analysis  of  distributed 

parameter  systems  is  also  found  in  the  areas  of  aerospace 

32 

vehicles  and  continuous  furnances  ,  panel  flutter  in 

p  7  28 

supersonic  f  1  o  w  ^  ,  vibrating  shafts  ,  and  electrical 

r  op 

power  transmission  lines0’  .  In  processes  where  temperature 

and  pressure,  say,  are  variables  it  is  obviously  necessary 
to  ensure  that  these  variables  are  maintained  within  given 
bounds  . 

Aerospace  re-entry  vehicles  use  ablative  shields 
to  protect  the  space  capsule  from  aerodynamic  heating. 

This  is  an  example  of  a  variable-domain  distributed  parameter 
system  in  which  the  domain  boundary  motion  depends  upon 
certain  system  variables  such  as  velocity,  evaluated  at  the 
boundary.  Continuous  furnaces  are  used  in  drying  and  heat 
treating  processes.  It  is  necessary  to  maintain  close 
control  of  the  temperature  distribution  of  a  uniform  material 
strip  inside  a  furnace  by  means  of  varying  the  velocity  of 

t 

the  transport  mechanism.  The  electrical  power  transmission 
line  is  an  example  of  a  fixed  domain  distributed  parameter 
system  where  it  is  required  to  closely  regulate  the  load 
voltage  and  generator  frequency. 


. 
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1.2  Different  Types  Of  Stability 

Basically,  there  are  two  main  types  of  stability 
which  are  covered  in  the  literature;  bounded  input-bounded 
output,  and  null  state  or  Liapunov  stability.  The  former 
type  has  been  investigated  by  several  authors.  Zames38 

used  functional  analysis  methods  to  solve  the  output 

stability  of  certain  nonlinear  time-varying  systems. 

16  2  2  25 

Others  *  have  used  the  maximum  principle  for  parabolic 

26  3 
equations  to  extend  the  work  of  Bellman*3  on  the  nonlinear 

13  1  R 

equation  of  heat  conduction.  Kastenburg  used  the 

method  of  comparison  theorems  to  solve  the  stability  problem 
for  a  class  of  nonlinear  feedback  control  systems  governed 
by  parabolic  partial  differential  equations. 

This  type  of  stability  will  not  be  discussed 
further  here.  In  this  thesis  only  Liapunov  stability  will 
be  investigated.  Suffice  it  to  say,  however,  that,  as 
shown  by  De  Figueiredo8,  bounded  input-bounded  output 
stability  can  be  deduced  under  certain  conditions,  from 
asymptotic  stability.  The  basic  concepts  involved  in 
Liapunov  stability  are  discussed  in  the  following  section. 

1.3  Liapunov  Stability 

a)  Lumped  systems 

For  a  system  of  ordinary  differential  equations, 

R  =  f (R)  , 

the  following  is  the  so-called  e,5  definition  of  stability 
in  the  sense  of  Liapunov; 


' 
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DEFINITION 


fl 


An  equilibrium  state  of  a  free  dynamic  system 

is  stable  if  for  every  real  number  e  >  0  there  exists  a 

real  number  S(e,t  )  >  0  such  that 

o 


i.  -  *e 


<  6 


<  e 


for  all  t  >  t 


imp! i e  s 

|| j> ( t ;R0 . t o )  -  Re 
where  R_0  is  the  initial  state  vector,  is  the-solution 
vector  of  the  above  system  of  differential  equations  and 
|| .  |  represents  a  norm  or  metric. 

If,  in  addition  to  the  above  definition,  the 
following  holds: 

||  $(t  ;R0  ,t0 )  -  Re  ||  0  as  t  +  “ 

then  the  equilibrium  state  Re  is  said  to  be  asymptotically 
s tabl e . 


The  idea  of  the  Liapunov  technique  is  to 
mathematically  describe  a  kind  of  distribution  of  stored 
energy  of  a  system.  It  is  the  purpose  of  the  direct 
method  of  Liapunov  to  determine  whether  this  energy 
decreases  toward  a  minimum  value  or  toward  zero.  If  this 
can  be  shown,  the  steady  state  of  the  system  is  stable 
or  asymptotically  stable,  respectively.  The  steady  state 
of  a  free  dynamic  system  described  by  ordinary  differential 
equations  is  defined  by  a  point  in  the  phase  diagram  of  the 
state  vari abl es . 

The  method  consists  of  choosing  a  scalar  function 
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V(jO  (Liapunov  function)  of  the  state  variables  which  has  a 
definite  sign.  If  the  sign  of  the  time  derivative  of  this 
function  is  semi -def i ni te  or  definite  and  of  the  opposite 
sign  then  the  steady  state  of  the  system  is  stable  or 
asymptotically  stable,  respectively.  Thus,  knowledge  of 
the  system  stability  is  found  without  having  to  solve  the 
transient  equations. 

The  disadvantage  of  this  technique  is  that  it 
yields  only  a  sufficient  condition  for  stability.  If  the 
sign  definiteness  criteria  are  not  met,  nothing  can  be 
concluded  about  the  stability  of  the  system.  When  this 
happens,  another  function  must  be  chosen  as  a  Liapunov 
function  and  try  the  technique  again. 

For  linear  stationary  systems  it  has  been  found 
that  for  global  asymptotic  stability,  the  following 
Liapunov  function  is  both  necessary  and  sufficient. 

V  (R_)  =  RTP  R 

where  IP  i  s  a  positive  definite  symmetric  matrix. 

Methods  for  determining  stability  of  these 
constant  linear  systems  are  well  established.  Many 
excellent  books  covering  the  subject  are  available, 
including  Hahn  *  * ,  Krasovski i ^  and  LaSalle  and  Lefschetz^0, 
to  mention  a  few.  Although  general  nonlinear  systems 
cannot  yet  be  treated  satisfactorily,  Liapunov  stability 
theory  has  helped  to  identify  various  time-varying  and 
nonlinear  models  with  well-established  qualitative 
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behavior.  Brockett,  in  his  survey  paper  discusses  many  of 

these  models  in  detail.  Many  of  the  design  techniques  of 

classical  control  theory,  such  as,  Nyquist's  criterion, 

1 8 

Root-Locus  and  Root-Contour  methods  ,  all  are  directly  or 

29 

indirectly  related  to  the  circle  criterion  ,  which  can  be 
derived  from  Liapunov  stability  theory, 
b)  Distributed  systems 

Stability  of  distributed  parameter  systems  is 
much  more  difficult  to  visualize  than  stability  for  lumped 
parameter  systems.  Instead  of  being  a  single  point  in 
state  space,  the  steady  state  now  consists  of  a  profile 
which  is  a  continuum  of  an  infinite  number  of  equilibrium 
points.  Thus,  the  idea  of  phase  planes  becomes  meaningless 
since  a  separate  phase  plane  diagram  would  be  needed  for 
each  of  these  equilibria. 

Now  stability  depends  on  the  proximity  of  the 
transient  solution  (that  is,  the  response  arising  after  a 
perturbation  about  the  steady  state  profile)  to  the  steady 
state  profile.  If  the  transient  response  remains  within  a 
finite  region  of  the  steady  state  the  system  is  stable 
about  that  equilibrium  profile.  Asymptotic  stability  means 

t 

that,  given  sufficent  time,  the  difference  between  the 
steady  state  profile  and  the  transient  profile  will  become 
as  small  as  desired.  Mathematically,  Liapunov  stability  of 
systems  described  by  partial  differential  equations,  such 


as 


=  f(U) 
9  t 


s 


' 
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is  defined  analogously  to  stability  of  lumped  systems.  Here 
U,  the  state  solution  vector,  is  defined  as  a  perturbation 
from  the  equilibrium  vector.  By  making  an  appropriate 
change  of  coordinates  the  equilibrium  vector  can  be  made  to 
coincide  with  the  origin  or  null  vector.  Thus,  the 
following  definition  is  sometimes  called  the  definition  for 
null  state  stability: 

DEFINITION 


An  equilibrium  state  U.e  =  0  is  stable  if  for 
every  real  number  e  >  0  there  exists  a  real  number  6(e,t0)  >0 
such  that 

|U(t0,X)  ||n  <  «(e.t„) 


implies  that 


U(t,X) 


£  e  for  all  t  >  t0 


If,  in  addition,  lu_( t ,?L)  0  as  t  ■+  °°,  the  null  state  is 

i  C 


n 


denotes  the  effect 


asymptotically  stable.  The  notation 
that  the  dependence  of  U_(t,X_)  on  the  position  X_  in  the 
spatial  domain  has  on  the  norm.  For  example,  the  norm  of 
the  state  vector  could  be  taken  as  the  L2  norm, 


lHU.Dln  =  j-l/n  dn  1 


The  most  general  extensions  of  Liapunov's  work 

39 

to  infinite  dimensional  systems  has  been  done  by  Zubov 
He  extended  the  previous  notions  concerning  a  Liapunov 


function  for  finite  dimensional  systems  into  the  more  general 
concept  of  a  Liapunov  functional  which  can  be  used  in  the 
infinite  dimenstional  space.  A  functional  is  a  correspondence 
which  assigns  a  definite  number  to  each  function  belonging  to 
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. 

. 
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a  certain  class.  Zubov's  main  results  are  summarized  in  a 

32 

theorem  in  Wang 

The  Liapunov  stability  theorem  for  infinite 

dimensional  systems  is  a  special  case  of  Zubov's  theorem. 

It  is  directly  analogous  to  the  Liapunov  theorem  for  systems 

described  by  ordinary  differential  equations11,  that  is, 

12 

THEOREM 

The  sufficient  condition  for  a  distributed 
parameter  system  to  have  an  asymptotically  stable  equilibrium 
state  is  that  there  exists  a  positive-definite  functional  V 
which  has  a  negative-definite  time  derivative  V  along 
solutions  to  the  system  equations.  Such  a  functional  is 
called  a  Liapunov  functional. 

Just  as  in  the  simple  case  of  the  n-di mens i onal 
problem,  the  main  problem  in  using  the  Liapunov  technique 
in  distributed  parameter  systems  is  the  determination  of 
the  Liapunov  functional.  Zubov  has  suggested  for  a  system 
of  strongly  parabolic  partial  differential  equations  that 
the  Liapunov  functional  have  the  following  form: 

V  =  uT(t,x)  u(t.x)  da 

It  is  obvious  that  this  is  the  distributed 
analog  of  the  previous  quadratic  form  where  the  matrix  P_ 
has  been  reduced  to  the  identity  matrix  I_.  The  fact  that 
using  a  weighting  matrix  £  in  the  Liapunov  functional  yields 
sharper  stability  results  will  be  shown  in  later  chapters. 


'■ 
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1.4  Methods  of  Determining  Liapunov  Stability  of 

Distributed  Parameter  Systems 


9 


In  the  analysis  of  systems  with  distributed 

parameters,  usually  the  system  is  first  approximated  by 

a  lumped  parameter  model.  Then  known  techniques  for 

9 

lumped  systems  are  used  .  Although  this  is  a  seemingly 

practical  approach,  it  does  not  always  lead  to  satisfactory 

results.  For  example,  in  the  design  of  feedback  controllers, 

the  control  law  may  be  stable  for  the  lumped  parameter 

approximation,  whereas,  the  steady  states  of  the  actual 

33 

system  may  be  unstable 

12 

Hsu  mentions  some  of  the  problems  involved 
with  approximation  with  reference  to  the  areas  of  nuclear 
reactors,  heat  exchangers  and  pneumatic  transmission  lines. 

In  general,  the  approximations  take  one  of  the  following 

-  35 

forms 

1)  Spatial  discretization:  This  results  in  a 
finite  system  of  continuous  time  ordinary 
differential  equations. 

2)  Time  discretization:  The  resulting  model 
consists  of  a  finite-dimensional  system  of 
spatially  continuous  ordinary  differential 
equations. 

3)  Space-time  discretization:  The  model  takes 
the  form  of  a  finite-dimensional  system  of 
difference  equations.  This  form  is  useful 


■ 


.. 


I 


'■ 

■ 


■ 


10 


for  digital  computation. 

4)  Spatial  harmonic  truncation:  In  cases 
where  the  state  functions,  over  a  finite 
spatial  domain,  have  a  band-limited 
spectrum  with  high-frequency  spatial 
harmonics  attenuated,  the  system  may  be 
approximated  by  a  finite-dimensional 
system  by  truncating  the  spatial  harmonics 
at  a  suitable  frequency. 

In  many  cases,  however,  the  inherent  distributed 

nature  of  a  problem  is  lost  by  making  approximating 

assumptions  about  the  behavior  of  the  system  variables.  For 

2 

example,  problems  involving  epidemics  .  It  is  obvious  that 
the  epidemics  are  distributed  over  a  spatial  domain  and, 
although  they  have  been  modelled  as  partial  differential 
equations,  in  some  instances  the  assumption  is  made  that 
the  epidemic  is  isolated  in  such  a  way  that  no  spatial 
dependence  arises.  Since  the  resulting  model  is  highly 
unrealistic,  its  use  could  easily  yield  misleading 
conclusions. 

To  avoid  problems  involved  with  the  various 
types  of  approximation,  it  is  desirable  for  studying  stability 
to  formulate  the  problem  directly  as  a  distributed 
mathematical  model^*^.  The  distributed  models  involve 
integral  and  partial  differential  equations  which,  in  general, 
are  very  difficult  to  solve.  The  best  single  technique 
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currently  available  for  tackling  these  problems  is  Liapunov's 
direct  method  whichhas  been  very  useful  in  qualitatively 
describing  their  behavior  without  having  to  solve  them10*33. 

1.5  Scope  of  the  Thesis 

In  this  thesis,  stability  of  distributed 
parameter  systems  is  studied  using  Liapunov's  second  or 
direct  method.  Two  techniques  are  presented  for  finding 
Liapunov  functionals  for  processes  which  are  described  by 
systems  of  nonlinear  partial  differential  equations.  The 
first  technique  involves  an  extension  to  the  work  done  by 
Berger  and  Lapidus^  which  uses  a  matrix  quadratic  form  in 
the  integrand  of  the  Liapunov  functional.  This  method  was 
suggested  but  not  attempted  by  Berger  and  Lapidus.  Two 
applications  are  provided  to  illustrate  the  technique 
developed.  The  technique  is  first  applied  to  the  catalyst 
particle  problem  and  then  to  the  tubular  reactor  problem. 

These  two  problems  are  characterized  by  the  possibility  of 
having  multiple  steady  states. 

The  second  technique  developed  is  completely 
analogous  to  the  well-known  variable  gradient  method.  This 
method  has  found  wide  use  in  studying  stability  of  systems 
of  ordinary  differential  equations.  An  algorithm  is 
developed  which  extends  the  theory  to  the  distributed  parameter 
case.  For  the  purpose  of  comparison,  the  same  examples 
treated  by  the  other  technique  are  used  to  illustrate  the 
variable  gradient  method. 


. 
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CHAPTER  2 


LIAPONOV  STABILITY  USING  A 
FUNCTIONAL  OF  QUADRATIC  FORM 


12 


2.1  Introduction 

In  this  chapter  a  technique  is  developed  for 

determining  stability  of  distributed  parameter  systems 

which  uses  a  Liapunov  functional  having  a  quadratic  form 

21 

as  the  integrand  .  In  general,  the  problems  which  can  be 
treated  by  this  method  have  system  equations  of  the 
following  form: 


3U(t,X)  =  f[U(t,X)l 

Ft 


(2.1) 


where  ]J(t,)0  is  the  perturbation  vector  and  £,  in  general, 
represents  a  second-order  vector  nonlinear  differential 
operator  over  the  spatial  domain.  As  mentioned,  the 
Liapunov  functional  has  the  form: 


(2.2) 


where  P_  is  a  positive-definite,  symmetric  matrix. 

This  method  is  an  extension  of  the  work  of 


4 

Berger  and  Lapidus  .  Although  they  suggested  the  use  of 


the  non-identity  matrix  P_,  they  never  attempted  it.  The 
development  of  the  method  will  be  illustrated  through  two 
examples  since  the  details  will  be  made  clearer  than  by 
speaking  in  general  terms.  The  two  examples  are  the  same 
as  those  treated  by  Berger  and  Lapidus  and,  therefore,  the 


.1 


. 


' 


analysis  roughly  parallels  theirs. 
2.2  Catalyst  Particle  Problem 
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Multiple  steady  s ta tes some t i mes  result  when  a 
first-order  chemical  reaction  occurs  inside  a  porous 
catalyst  particle  in  the  presense  of  heat  and  mass 
diffusion.  A  physical  interpretation  of  these  multiple 
equilibria  may  be  obtained  from  the  following  diagram 

37  w 

from  Weisz  and  Hicks  .  The  curved  line  in  the  sketch 


Figure  1 

illustrates  the  functional  relationship  between  the  rate  of 
heat  production  with  temperature.  The  straight  line, 
representing  the  rate  of  heat  removal  by  conduction  for  a 
given  temperature  differential  between  adjoining  volume 
elements,  intersects  the  curve  up  to  three  times.  Thus, 
equating  heat  removal  and  production  rates  results  in  up  to 
three  steady  state  solutions. 

The  point  b  is  metastable.  An  initial  temperature 


. 


. 
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condition  below  b  results  in  stabilization  at  point  a,  and 

initial  conditions  corresponding  to  a  point  above  b  sends 

the  system  to  condition  c.  Thus,  in  observing  reaction 

rates  as  a  function  of  rising  temperature,  a  discontinuous 

temperature  change  may  obtain  when  a  metastable  condition 
// 

is  passed 
MATH  MODEL 

For  the  purpose  of  analysis,  the  catalyst  particle 
is  assumed  to  be  in  the  form  of  a  slab.  The  following 
dimensionless  mass  and  heat  balance  equations  are  associated 
with  the  first-order  reaction  within  the  particle. 


£X  =  llx  -  4>2y  exp 
at  a  x2 


y(z-i) 


N  3^z  ^  3 2 z  +  3 d> 2 y  exp 

3 1  ax7 


Y(z-l) 


-] 


with  the  boundary  conditions  given  by 


3y(  t  ,0)  =  3z_[t,0)  =  0,  y  ( t ,  1 )  =  z  ( t ,  1 )  =  1 
3x  3x 


(2.3) 

(2.4) 


(2.5) 


The  usual  symbols  are  used  here  and  their  meanings  are 
stated  in  the  List  of  Symbols  on  page  iv.  The  detailed 
derivation  of  these  equations  can  be  found  in  APPENDIX  2. 

At  steady  state,  equations  (2.3)  and  (2.4)  hold 
but  with  the  time  derivatives  becoming  zero.  Denoting  the 
solutions  of  the  resulting  steady  state  equations  by  y*  and 
z*,  an  adiabatic  balance  at  this  equilibrium  condition  yields 


■ 
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. 
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z*  =  1  +  3 ( 1-y*)  (2.6) 

To  analyse  the  stability  about  the  equilibrium  states, 
define  the  following  set  of  perturbation  equations  by 
letting 

y  =  y*  +  u 
z  =  z*  +  v 


(2.7) 


and,  writing  the  nonlinear  reaction  rate  term  as 


f(y,z)  =  <f>2y  exp 


y(z-l) 


(2.8) 


then,  (2.3)  and  (2.4)  yield,  respectively 


3u 


at  ax 


t  -  f(y,z)  +  |4* 


ax 


(2.9) 


N  _av 

at 


a2v  .  or/..  ,  a2z* 


ax; 


+  3f(y,z)  + 


ax- 


(2.10) 


with  the  boundary  conditions  becoming 

Mt.O)  _  _  U(t,l)  =  v ( t ,  1 )  -  0 

ax  ax 


(2.11) 


Linearizing  the  nonlinear  reaction  rate  term  about  the 
steady  state  profile  with 

v 

y=y* 

z=z* 

=  f (y * )  +  fyU  +  fzv  (2.12) 

and  using  the  steady  state  equations 


f(y»z)  =  f(y*,z*)  + 


9f (y »z) 

ay 


u 

y=y* 

z=z* 


3f(y,z) 

az 


. 


■ 


. 

. 
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||T  -  f(y*)  (2.13) 

|p-*  -  -$f(y*)  (2.14) 

yields  the  following  linearized  perturbation  equations 

-  v  -  V  (2-15) 

II  -m|p  +  MBfyU  +  MBfzv  (2.16) 

where,  for  convenience,  M  denotes  1/N. 

Guided  by  the  suggestions  of  Zubov  for  strongly 
parabolic  equations,  Berger  and  Lapidus  use  the  following 
Liapunov  functional 


V  =  X  f  UT(t,X)  U(t,X)  dfi 

2  J  0  ~  -  -  — 

(2.17) 

where 

^  U  If 

T 

(2.18) 

u  -  [u  V ]  , 

the  state 

vector.  In  this  case  (2.17)  becomes 

V  =  i  f  (u2  +  v2 )  dx 

2  J  U 

(2.19) 

As  pointed  out  at  the  outset,  the  purpose  of  this 
chapter  is  to  investigate  the  use  of  a  non-identity  matrix 
in  the  Liapunov  functional.  For  this  problem,  equation 
(2.2)  takes  the  following  form 

v  -  i  f  UT(t,x)  P  U(t,x)  dx 

2  J  A 


(2.20) 


. 
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£  is  a  positive-definite  symmetric  matrix,  written  for 
convenience  as 


Pi  P3 
Pa  P2 


(2.21) 


Expanding  (2.20)  with  (2.18)  and  (2.21)  yields 
V  =  JL  f  (p  u2  +  2p  uv  +  p  v2)  dx  (2.22) 

2  Jo  1  3  2 

Since  P_  is  positive-definite,  V  is  positive-definite. 

From  (2.22) 

v  u!r  +  ^(ufr  +  ^  +  v!r]dx  <2-23> 

o 

Substituting  (2.15)  and  (2.16),  (2.23)  becomes 


p  u 


3  2  U  „  ..  2 


1  3x2 
3  2  v 


-  PJUzfy  -  PjUVfj, 


+  p  Mu-2 — 1  +  p  MBu  f.;  +  p  MBuvf7 
‘  3x2  3  y  3  z 

s2  ii 

+  p  v- 

‘  3x‘  3  3 

2*  dx  (2.24) 


3  3  x 2  3 

'•5-ii  -  p  uvfv  -  p  v 2 f 7 
a  ax2  a  y  Ha  z 

3  2  v 


+  p  M v- — -  +  p  MBuvfv  +  p  MBvzf7 
2  3x2  2  y  2  Z 

Using  integration  by  parts  on  second-order  terms,  (2.24) 
becomes 


-  P  (— )  (— )  -  P  uvf  -  p  V2f 

3  3  X  3  X  3  y  a  L 

-  p  M (— ) 2  +  p  M$uvf  +  p  Mg  v2  f. 

2  3X  2  y  2  ' 


dx  (2.25) 


■ 
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Consider  the  following  inequality 

0  (m  ±  n)2  =  m2  ±  2mn  +  n2 
Thus  ±2mn  _<  m2  +  n2 

or,  in  general 

±amn  _<  | a | | mn |  £  |aj (m2  +  n2) 

2 

Using  (2.26)  in  (2.25)  and  rearranging  yields 

V  £  Pu2(p  M6fy  -  P^y  +  “|  P  3M3f  z  “  Pxfz  -  P^y 
J  o 

+  J  V( P2MSfz  -  Pjf2  +  ||P3MBfz  -  P3fz  -  P3fy 
+  /o'(^)J(-P,  +7lP3(H+l)|)  dx 
+  f  ’( fx)J(~P2M  +  Tlp3(M+1)  I  )  dx 
To  simplify  this  expression  let 

Fu  =  ( p 3 H B  -  p3)fy  +  yI(P3MB  -  P3)fz  +  (P2MB  -  Ps 
Fv  =  (p  MB  -  P  )f2  +  7I  (p  MB  -  P  )fz  +  (p  MB  -  p 
Then  (2.27)  becomes 

V  £  J  1  u 2 F u  dx  +  J  v  2  F  y  dx 
-  /.‘(^‘(P,  -  f<H+l>IP3l)  dx 

fj)2(pi"  '  t<m+1>Ip3I>  dx 

g 

Now  consider  an  inequality  credited  to  Gavalas  . 


(2.26) 

+  p^MBfy | ) dx 
+  p  MBfy | ) dx 

(2.27) 

)fyi 

(2.28) 
)fyl 

(2.29) 


(2.30) 

From  the 


- 
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i den  ti ty 


u(t,x)  =  -  f  iy.(t »x)dx  »  for  u ( t ,  1 ) =0 

Jx  3x 


it  follows  that 

u2(t,x)  = 


r 


3  u  ( t ,  x ) 


3x 


dx 


and  with  the  Schwarz  inequality 


31 


or 


u2(t,x)  <  f  (|^)2 
J  x  3x 

u2(t,x)  <  (i-x)  f  (iy.) 

J  x  3x 


dx  f  ( 1 ) 2  dx 
J  x 


dx 


(2.31) 


Similarly 


v 2 ( t * x )  <  d-x)  rw  dx 

J  x  3x 


(2.32) 


With  (2.31)  and  (2.32)  equation  (2.30)  can  be  written 

dx 


V  <  I  F 


u 


3  u  ■  2 


a- 

( l-x)  f  (|^-)  dx  I  dx 


3x 


J 


i  (ll)2!p>  '  4-<m+1> !p3']  dx 

J'(f-V  P,"  -  i(M+l)|p,|]  dx  (2.33) 


rearranging  finally  yields 


V  <  1  I  (£“)*  dx !  ; - { p  -  — ( M+ 1 ) | p  |)  +  /  (l-x)F  dx  ! 
p  ix  M  1  2  3  0  u  j 


'  \  I 

3  u  \  2  A  v  \ 


1  \i  r 1  \ 

(1^)  dx;  -{p  H  -  ~( M+ 1 )  |  p  |)  +,  (l-x)F  dx 
3x  j\  2  2  3  0  v  J 


-  0 


(2.34) 


' 


. 


. 
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It  is  apparent  from  (2.34)  that  two  criteria 
must  hold  simultaneously  to  give  a  sufficient  condition 


asymp to ti c 

stability. 

These  criteria  are 

(a) 

(VxjFu 

J0 

dx  <  p(-  ~( M+l ) | p s | 

(2.35) 

(b) 

f\  l-x)Fv 
-^0 

dx  <  P2 M -  j(H+l)|p3| 

(2.36) 

DISCUSSION  AND  RESULTS 

The  derivation  here  differs  slightly  from  the 
work  of  Berger  and  Lapidus.  They  separated  the  integral  of 
Fu  into  two  sets  corresponding  to  intervals  where  Fu>0  and 
intervals  where  Fu<0.  The  integral  of  Fv  was  not  similarly 
manipulated,  since,  in  their  case  P_=l,  the  identity  matrix, 
and  thus,  Fv>0  for  O^x^l .  For  their  special  case  of  N=l, 
in  place  of  (2.30)  they  obtained 

V  <  -  f  [(ii)2  +  (— )2]  dx  +  f  v2  F  dx 
~  J0  3x  8x  Jo  v 

3" 

I F u | u 2  dx  - 

it 

xe[  0  ,1] 

In  the  next  step  Berger  and  Lapidus  throw  out  the  last  term 

i 

for  no  apparent  reason  except  that  it  strengthens  the 
inequal i ty . 

As  will  be  shown  the  retention  of  all  terms  in 
(2.30)  is  advisable  when  a  weighting  matrix  is  used  in  the 
Liapunov  functional.  It  was  found  that  both  Fu  and  Fv  could 


Fulu 

'xd  [  a"  ,  £"] 


dx 


(2.37) 


- 


. 


-■ 
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oscillate  from  positive  to  negative  with  x  and  that 
restricting  x  to  only  those  intervals  where  the  functions 
were  positive  resulted  in  cases  from  which  no  conclusion 
about  stability  could  be  drawn.  However,  for  the  same  cases 
retention  of  all  terms,  that  is,  not  restricting  x,  allowed 
conclusions  of  asympotic  stability  to  be  made. 

The  sufficient  conditions  for  asymptotic  stability 
which  result  from  re strict ing  x  to  only  those  intervals 
where  Fu>0  and  Fy>0  are 


( a  )'J 

( 1-x) F  dx  <  p  -  — ( M+l ) | p  | 

Fu>0  U  12  3 

(2.38) 

(b)j 

(l-x)F  dx  <  p  M  -  — ( M+l ) | p  | 
Fv>0  v  2  2  3 

(2.39) 

It  is  obvious  from  the  two  sets  of  suff i ci ent  conditions 
(2.35),  (2.36)  and  (2.38),  (2.39)  that  the  first  set  would 
be  the  easier  to  satisfy.  For  example,  split  up  the  left 
hand  term  in  (2.35)  and  compare  it  to  the  left  hand  term  in 
(2.38).  That  is,  compare 


with  just 


(l-x)Fu  dx 


=f  (l-x)Fu  dx  -f  ( 1  -  x )  |  Fu  |  dx 
JFU>0  JFU<0 


/  (l-x)F 
JF.,>0  U 


dx 


Thus,  in  trying  to  satisfy  the  respective 
sufficient  conditions  for  asymptotic  stability  it  will  be 
more  difficult  in  the  case  of  condition  (a)’  since  the  left 


' 


. 


p.  I 

■ 


' 

■' 
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hand  side  has  lost  a  term  which  tends  to  reduce  its 
magnitude  while  the  right  sides  remain  unchanged. 

As  mentioned  previously,  this  system  has  three 
possible  steady  states,  two  of  which  should  be  stable  and 
one  metastable.  By  using  the  analysis  outlined  above  for 
the  special  case  of  N=l,  Berger  and  Lapidus  were  able  to 
conclude  only  that  the  first  steady  state  is  asymptotically 
stable.  By  using  a  weighting  matrix  and  by  retaining  all 
terms  of  inequality  (2.30)  both  steady  states  one  and 
three  were  shown  to  be  stable  for  various  values  of  N.  Some 
of  the  results  are  shown  in  table  1  along  with  the  results 
of  Berger  and  Lapidus.  In  all  values  of  N,  one  result  is 
shown  where  dropping  the  above-mentioned  terms  yields  no 
conclusion  about  stability. 

The  IBM/APL  360/67  computer  system  with  direct 
access  time  sharing  facilities  was  extremely  efficient  for 
handling  the  computations  required  for  this  problem. 
Although  suitable  elements  of  the  matrix  £  can  be  found 
only  by  trial  and  error  methods,  the  APL  system  made  this 
task  s trai ght-forward . 

As  the  identity  matrix  enabled  stability 

i 

conclusions  to  be  made  for  steady  state  one  for  most 
values  of  N,  it  was  always  used  as  an  initial  trial  for 
each  steady  state.  Then,  if  that  failed  to  produce  a 
conclusion  about  stability,  the  following  general  procedure 
was  implemented. 


■■ 


' 
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These  values  were  obtained  using  conditions  (a)’  and  (b)\ 
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Without  loss  of  generality,  let  p  =1.  The  values  of  p  and 

1  2 

P3  are  then  chosen  arbitrarily,  keeping  in  mind  that  must 
be  posi  tive- def  i  ni  te  .  Inequalities  (2.35)  and  (2.36)  are 
tested.  If  they  are  not  satisfied,  the  value  of  p  (or  p  ) 

2  3 

is  changed  and  the  tests  are  repeated.  If  the  tests  fail 
again  one  can  compare  the  amounts  by  which  the  inequalities 
failed  to  be  satisfied  and  use  these  amounts  to  immediately 
determine  whether  the  direction  of  the  change  chosen  for  P2 
(or  p^)  was  correct.  That  is,  if  the  inequality  held  in 
reverse,  the  amount  by  which  the  left  hand  side  exceeded 
the  right  hand  side  serves  as  a  measure  of  convergence  or 
lack  of  convergence  toward  a  satisfactory  P_  matrix.  It 
was  sometimes  necessary  to  also  vary  the  third  element  to 
finally  arrive  at  such  a  matrix. 

2.3  Tubular  Reactor  Problem 

In  this  section  the  method  developed  earlier  is 
further  illustrated  by  means  of  a  second  example.  A  simple 
first-order  chemical  reaction  which  takes  place  in  a 
tubular  reactor  is  analysed  to  determine  the  stability  of 
the  steady  states.  The  fluid  in  the  reactor  is  assumed  to 
be  in  turbulent  flow  such  that  radial  effects  may  be 
neglected.  Thus,  the  only  diffusion  which  occurs  is  in  the 
axial  direction. 

MATH  MODEL 

The  normalized  mass  and  heat  balance  equations 


are 


' 

i 


■ 

> 


. 


25 


(2.40) 


(2.41) 


wi  th 


(2.42) 


where,  again,  the  symbols  are  as  defined  in  the  List  of 
Symbols  on  page  iv. 


As  in  section  2.2,  generalized  perturbation  is 


used  here  to  determine  the  stability  of  the  steady  states. 

These  generalized  results  are  later  reduced  to  the  special 

case  of  the  adiabatic  perturbation  and  the  results  so 

4 

obtained  are  compared  to  those  in  . 

Consider  the  nonlinear  term 


(2.43) 


f(y,z)  =  y  exp ( -Q/ z ) 


Linearizing  this  about  the  steady  state  y*,z*  yields 


f(y,z)  =  f(y*»z*)  +  gy  +  —  (?-?*) 


(2.44) 


At  steady  state,  (2.40)  and  (2.41)  become 


(2.45) 


ilz*  =  -  vf (y* , z*) 

3s  2  ^  3s 


(2.46) 


;  .  j  t-n-n 


~  ■ 

' 


26 


So,  if  the  perturbations  in  y  and  z  are  denoted  as  before  by 
u  and  v,  respectively,  the  transient  equations  can  be  written 
as 


3_u  _  32u 

at  "  Ts* 

3v  _  3 2  V 
3 1  =  Is7 


N  — 
P  3  S 


32y* 

37r 


N 


3v_ 

P3s 


+ 


3  2  z* 

dU 


3  V* 

Np—  +  af(y*,2*)  +  afyu  +  afzv 

(2.47) 

3  z  * 

Npg"s  +  vf(y*,z*)  +  vfyu  +  vfzv 

(2.48) 


combining  the  last  four  equations  yields 


3u 

3  2  u 

-  N  — 
P3s 

3 1 

3s* 

3  v 

3  2  v 

-  N  — 
P3s 

3 1 

3sJ 

+  afyU  +  afzv 


+  vfyU  +  vfzv 


(2.49) 


(2.50) 


with  the  boundary  conditions 

u(t,0)  =  v(t,0)  =  0  (assumed) 

and 


(2.51) 


=  =  0  (from  <2-42)) 


To  make  analysis  easier,  the  following  transformat ions 
will  be  used 

r n  si 

(2.52) 


u ( t , s )  =  £(t,s)  exp 
v ( t , s )  =  n(t,s)  exp 

L  2  J 

These  relations  transform  (2.49)  and  (2.50)  into  the 
following 

£  '  : &  +  -  "h  +  onf. 


(2.53) 


3t  3s 
,  2 


(2.54) 


-  |_4  +  n(vf. 


N2 


3s 


41  +  v5fy 


9_n 
3 1 


(2.55) 


. 


. 


" 


•  S'  I  ■ ,  ' 
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with  the  boundary  conditions 


£ ( t ,  0 )  =  n(t  ,0)  =  0 


N, 


9s ( t , 1 )  =  -  2  ^ ( t  *  1 ) 


3n 


N, 


9  s  (  t  ,  1  )  =  — ^T|  (  t  ,  1  ) 


(2.56) 


As  i n  section  2.2,  let 


V  =  —  |  UT(t,s)  P  U(t,s)  ds 
2  -'0 


(2.57) 


where 


and,  as  before 


U  = 


P  = 


S(t,s) 
n( t  ,s ) 

Pi  P  3 
P3  P2 


Differentiating  (2.57)  with  respect  to  time  yields 


V  =  f[ P.sff  ♦  P3U|t  +  nff)  +  P2nf^]ds  (2.58) 

-'O 

Substituting  (2.54)  and  (2.55)  into  (2.58)  yields 


flf  32  r  N* 

v  =  i  p.5{a^  '  7  «  +  a5fy  +  anV 

,2 


+  P 


5{gn  .  «pn  +  +  vny 


n( 


2  r  N 


+  P  n( 


ill 

8S  2 

32n 


_  _P£  +  a£f  +  apf  } 


N 


—  n  +  v£f  +  vnf  } 
2  9s2  4  y  1  zJ 


ds  (2.59) 


. 
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Integration  by  parts  of  some  of  the  terms  in  (2.59)  yields 
the  following 


3s 


ds 


=  € 


i£ 

3s 


ds 


■  -  f ' f>' 

■  f  (|f)=  I* 


ds 


(2.60) 


Simi 1 arly 


also 


N 

=  -JV(t,l) 
2 


(2.61) 


=  -^5(t.l)n(t,l)  -  fQ  <ff)  ds 

(2.62) 


Equation  (2.59)  can  now  be  written 

N 
2 


V  =  -  P  ^52(t,l)  -  P.N  S(t,l)n(t,l)  -  p  -En5(t,l) 

19  3  p  *2 


p  [’(!£) 2  ds  -  p J\^)2  ds  -  2P,  r(|£)(f)  ds 

iJ0  3s  2^o  3s  3^o  3s  3s 

ft'z2[pi(aTy  -  ^  +  p3vfylds 

/  1,2  [P2(vfZ  '  +  P3“f2  ]  dS 


N, 


J  Cn[piafz  +  P3(vfz  -  +  afy)  +  p2vfy]  ds  ^2*63^ 


Consider  the  following  equality 


£(t,s)  =  fS(^-)  ds 

Jo  3s 


■ 


. 
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therefore 


S2(t,s)  = 


(^)  ds 


9  S 


(2.64) 


Now,  with  Schwarz's  inequality 


? 2 ( t  ,s ) 


S  2 

(1)  ds 


\ 


V  ,  M 

<  s  /  (2£)  ds 

"  9s 


fsm 

9s 


ds 


(2.65) 


Similarly 


n2(t,s)  £  s  I  (!H)  ds 
Jo  8s 


(2.66) 


The  product  term  can  be  handled  using  the  inequality  used 
in  section  2.2,  that  is 


±a£n  <  | a | | |  <  I  a  |  U2  +  n2 ) 


(2.67) 


Upon  substitution  of  ( 2 . 56 ) , ( 2 . 66 )  and  (2.67)  into 
(2.63)  yields 

V  <  -  pi^2(t,l)  -  p2-^n2(t,l)  -  P3Np£(t,l)n(t,l) 


j  (|f)2  ds 


N 


-  P  +  s  p 


(afv  -  -e0  +  P3vfy  ds 


i  ^  y 


2 1 P3 1  +  l  fiP.afz  +  P3(vfz  -  +  “V  +  p2vf: 


|  ds 


I  (f11)2  ds 
Jo  9s 


-  p2  +  I  s  p 


N 


(vf7  -  -E-)  +  p 3 ctf z  ds 


2 '  v  1  z 


+  2  I P  3 1  +  l  f 1 P iafz  +  P3(vfz  -  *  afy)  +  P2vfylds 

(2.68) 


To  simplify  this  expres s i on, 1 e t 


- 


' 


30 


Np 

p,(afy  •  T>  +  P=vfy 


and 


+  j|p  ofz  +  p  (vf2 


N, 


+  <*fy)  +  P  2 vf y |  (2.69) 


Fn  =  P  .(vf. 


N, 


■)  +  P  af_ 

3  Z 


+  tI  P  af 7  +  P  • 


ND 

—  +  afv)  +  P  vfv|  (2.70) 

2  J'  i  j 


Then 


V  <  -  -  P2^V(t,l)  -  p  H  E(t,l)n(t,l) 


(  <f7>’  d! 


*  '  'U1'  "■ 


-  p.+  zlp,l  +  J0  sFe  ds 

-  P2+  2 1 P 3 1  +  l  sFn  ds 


(2.71) 


Consider  only  the  first  three  terms  of  (2.71) 


-£p  C2(t,l)  -  N  p  £( t ,  1 )n ( t ,  1 )  -  p  — n2(t,l) 

2  i  p  3  z  2 


N, 


‘Pi 


£2 ( t , 1 )  +  2^-C(  t  ,l)n(t,l)  +  i_n2(t,l) 

Pi  Pi 


(2.72) 


Since  matrix  £  is  positive-definite, 


p  p  -  p2  >  0  and  p  >  0 

K 1  K2  K  3  K 1 


Hence 


p  <  +/  p  p 
3  1  2 


P  <  +/p  p 
3  1  2 


(2.73) 


By  (2.73), 


the  terms  in  (2.72)  can  be  written 


- 
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N„  N„ 

-  JPi£2(t,l)  -  NpP35(t,l)ti(t,l)  -  P2~V(t,l) 


NP 

<  2  P, 


/j)  2 

5(t,l)  -  -A-n(t.l) 

/P, 


(2.74) 


And,  since  Np>0,  it  is  clear  that  the  first  three  terms  of 
(2.71)  are  negative-definite.  Hence,  in  order  for  V<0,  the 
following  conditions  are  sufficient 


f  sF  ds  <  p  -  2 | p  | 
J  0  C  1  3 


/  sF  ds<p  -  2 1  p 
'o  n  2  1  3 


(2.75) 

(2.76) 


For  the  special  case  of  adiabatic  perturbation, 
the  following  relation  holds 


3S  -  -n 


(2.77) 


When  equation  (2.77)  is  used  in  (2.71),  the  following 
sufficient  condition  for  asymptotic  stability  is  obtained 

fo  +  Frq)ds  <  Jz  +  P2  "  2  I  P 3  I  (  “-™2----)  (2.78) 

DISCUSSION  AND  RESULTS 

In  solving  for  the  steady  state  profiles  the 
following  parameter  values  were  used 

3  =  0.40,  N  =  30,  Q  =  30  and  a  =  2.25xl013 

4 

which  are  the  same  values  that  were  used  in  .  The  profiles 
were  generated  on  an  APL/360  direct-access  system*  using  a 

★ 

For  the  specific  functions  used,  see  the  APPENDIX 


\ 


* 


. 


; 


. 


' 
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double  precision  variable  step-size  Runge-Kutta  method. 

The  integration  was  carried  out  backwards  toward  the  inlet 
of  the  reactor  until  the  so-called  Danckwert  boundary 
condition  was  satisfied  there.  The  normalized  output 
temperatures  which  satisfied  these  input  conditions  are 
shown  below 

Steady  state  no.  Normalized  output  temperature 

1  1.053 

2  1.3166  . 

3  1.399999 

The  same  algorithm  which  was  roughly  outlined  in 

the  previous  section  was  used  here  to  find  the  elements  of 

the  matrix  £  for  both  steady  states  2  and  3.  Steady  state  1 

proved  to  be  asymptotically  stable  by  using  the  identity 

matrix.  For  steady  state  2  the  method  failed  to  converge 

to  a  matrix  £  which  would  satisfy  the  sufficient  conditions 

(2.75)  and  (2.76).  Although  this  is  not  a  confirmation  of 

37 

the  fact  that  steady  state  2  is  unstable  ,  it  does  show 

that,  in  this  case,  the  method  does  not  yield  matrices 

which  would  indicate  that  the  steady  state  is  stable.  The 

results  are  contained  in  Table  2. 

The  results  for  the  adiabatic  tubular  reactor, 

obtained  using  the  sufficient  condition  (2.78),  are  shown 

4 

in  Table  3.  A  comparison  of  these  results  with  those  in 
shows  that  the  use  of  a  non-identity  matrix  £  in  the 
Liapunov  functional  has  been  successful  in  proving  that 
steady  states  1  and  3  are  both  asymptotically  stable, 
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whereas,  the  use  of  an  identity  matrix  proved  only  that 
steady  state  1  is  stable.  Table  2  shows  that  the  tubular 
reactor  subjected  to  generalized  perturbations  in 
concentration  and  temperature  also  exhibits  two 
asymptotically  stable  steady  states. 


TABLE  3 


STEADY 
STATE  # 

p 

p, 

MATRIX 

P,  P, 

2  3 

INEQUALITY  (2.78) 

COMMENT 

L.H.S. 

R.H.S. 

1 

1 

1 

0 

-513.52 

7.25 

Stab! e 

2 

1 

1 

0 

1444.6 

7.25 

No  conclusion 

3 

1 

1 

0 

7346.5 

7.25 

No  conclusion 

3 

1 

8 

-3.5 

-757.37 

-36.45 

Stab! e 

I 
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CHAPTER  3 

VARIABLE  GRADIENT  METHOD  FOR 
CERTAIN  DISTRIBUTED  PARAMETER  SYSTEMS 

3.1  I ntroduc ti on 

The  objective  of  this  chapter  is  to  present  a 
method  analogous  to  the  variable  gradient  method  for  lumped 

9 

systems  proposed  by  Schultz  and  Gibson  .  This  method  has 
been  applied  to  a  few  distributed  parameter  systems.  Further 
work  may  yield  successful  results  for  other  systems  as  well. 

3.2  Variable  Gradient  Method  For  Lumped  Systems 

In  this  section,  the  variable  gradient  method  as 
applied  to  systems  described  by  ordinary  differential 
equations  will  be  briefly  reviewed.  The  purpose  of  this 
review  is  two-fold.  One,  to  make  the  thesis  somewhat 
self-contained  and,  two,  to  set  the  tone  for  the  extension 
of  the  method  to  distributed  parameter  systems  which  follows. 

The  main  difficulty  in  applying  Liapunov's  direct 
method  has  always  been  in  finding  Liapunov  functions  in  a 
straightforward  manner.  The  variable  gradient  method  is 
based  on  the  fact  that  if  there  exists  a  Liapunov  function 
capable  of  proving  asymptotic  s tabi 1 i ty  of  a  parti cul ar  system, 
then,  the  gradient  of  this  function  also  exists.  The 
introduction  of  a  variable  gradient  somewhat  reduces  the 
emphasis  on  the  ingenuity  and  experience  of  the  investigator. 
This  will  be  clear  as  the  discussion  proceeds. 


1 


" 


- 

■ 


' 


. 


36 


with 

Assume 

then 

or 

where 

From  VV 

I f  the 

hoi ds  , 
Choose 


Consider  the  system  of  equations 
R  =  f(R) 


(3.1) 


f  (0)  =  0 

a  tentative  Liapunov  function 
V  =  V (R ) 


(3.2) 


•  av  .  x  9V  . 

V  =  —  r  +  —  r  + 

9r,  i  9r  i 

V  =  (VV)TR  =  VV-R 


VV 

av 

9  V 

—  »  *  •  • 

9V 

§ 

.  9r  i 

9r  2 

9r 

,  find  V  as 

-  [vv.> 
fol 1 ows 

VV  ,  ... 

2 

,  vv 

9  V  . 

+  -  r 

3rn  " 


iT 


n 


(3.3) 


(3.4) 


r- 

=  VV 
JQ  — 


dR 


(3.5) 


stipulation  that 

9 ( VV  . )  9 ( VV  . )  .  .  .  0  „ 

— l  =  — ^  j  ;  =  l,2,...,n  (3.6) 


3rj 


3r. 


then  V(R^)  is  independent  of  the  path  of  integration 

the  simplest  path  as  follows 

V(R)  -I  1  VV  (u  ,0,...,0)  du 
“  J  o  1  1 

+r  2  W  (r  ,u  ,0,...,0)  du  +  ... 

J  „  2  12  2 


+  n  vv  (r  ,r 

n  i  2  in 


, u  )  du 


n 

(3.7) 


*,  U  -  f.K*  {t  •  (r  ' 
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Thus,  the  problem  of  determining  a  Liapunov  function  which 
satisfies  Liapunov's  theorem  has  been  transformed  into  the 
problem  of  determining  the  gradient  of  V  such  that  condition 
(3.6)  holds. 


The  following  procedure  is  used  to  evaluate  the 
Liapunov  function. 

1.  Assume  an  arbitrary  column  vector  for  VV  such  as 


VV  = 


a 

n  ri 

+ 

a  12  r  2 

+ 

...  + 

a  m  rn 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a 

iri 

+ 

a  2  2 

+ 

...  + 

annrn_ 

(3.8) 


The  a^'s,  in  general,  consist  of  a  constant  part  and  a 
variable  part  which  is  a  function  of  the  state  variables. 
Without  loss  of  generality,  a^  variable  parts  are 
functions  only  of  r^  except  for  ann  which  is  only  a 
cons tant( usually  2)  to  simplify  the  proof  that  V(RJ  =  C 

7 

(a  constant)  represents  closed  surfaces  . 

a 

2.  Determine  V  from  VV_  by  (3.4). 

3.  Constrain  V  to  be  at  least  negative  semi -def i ni te 

4.  Set  a  =  constant,  as  above,  and  determine  the 

nn 

remaining  a.,  by  using  relation  (3.6). 

5.  Recheck  V  to  be  sure  that  step  4  has  not  altered 
the  constraints  of  step  3. 

6.  Determine  V  from  equation  (3.7). 


' 


' 
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3.3  Extension  of  the  Variable  Gradient  Method 

to  Distributed  Parameter  Systems 


In  the  case  of  distributed  parameter  systems  the 
Liapunov  function  becomes  a  Liapunov  functional.  Now  the 
vector  system  equation  has  the  form 


=  f ( t , R )  (3.9) 

at 


since  R  =  R_(  t  ,x ) 

where  x  is  the  spatial  variable(in  the  one-dimensional  case). 
Now,  if  we  assume  the  existence  of  a  Liapunov  functional  V 
in  the  following  form 


(3.10) 


where  ft  represents  the  spatial  domain.  Then,  differentiating 
with  respect  to  time  yields 


V 


(3.11) 


assuming  that  8  is  independent 
be  written  as 

iH  -  M  iLi  +  M  lH2 

at  ’  ar:at  ar2at 


of  time.  The  integrand  may 


lnn 

arnat 


(3.12) 


or 


a_H 

at 


(vh)t^ 

—  at 


(3.13) 


Thus,  if  we  assume  an  arbitrary  form  for  VH  it  is  possible 
to  evaluate  M  and  H  in  a  method  similar  to  the  one 

at 

outlined  earlier  for  lumped  systems.  An  analogous 


■ 
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algorithm  to  use  to  evaluate  V  and  V  from  VH_  might  be  as 
fol 1 ows 

1.  Assume  VH  has  the  form 


VH  = 


auri  +  +  •••  +  s,nr 


12  1  2 


i  n 1  n 


alr1+a22+...+  annrn 


(3.14) 


2. 


—  =  ( VH ) T  f- 

at  at 


(3.15) 


3.  Constrain  V  =  /  —  dn 

to  be  at  least  negative  semi -def i ni te .  Do  not  merely 

constrain  —  <  0  since  this  would  be  much  too  restrictive. 
8t  = 

4.  As  before  set  a  to  some  constant  and  use  the 

nn 

following  n.(-n,~.l-).  curl  equations  to  determine  the 


remaining  a-jj 


— VHi )  =  ;  i ,j  =  1,2 . n  (3.16) 

3  r  -j  ari 


5.  Recheck  V  since  step  4  may  have  altered  the 
constraint  mentioned  in  step  3. 

6.  Determine  V  from 


V  -Jn  H  dfi 


(3.17) 


where 


R  j 

H  =  /  “(VH)  dR 
0 


Again,  the  simplest  path  of  integration  is  used,  that  is 


. 
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,0)  du  j 

. . .  ,0 )  du2  +  ... 

• »rn-i »un)  dun  (3.18) 


3.4  An  Example:  The  Catalyst  Particle  Problem 


The  method  outlined  above  is  first  applied  to  the 
catalyst  particle  problem  discussed  earlier.  Starting  from 
the  linearized  perturbation  equations  (2.15)  and  (2.16) 


—  =  ili  -  fvu  -  fzv  (3.19) 

3t  ITx7  y  2 

11  =  +  M$f  u  +  M$f 7 v  (3.20) 

a t  ax2  y  1 

Proceed  as  foil ows  : 

Assume 

VH  = 


a n  u  +  ai2  v 


a21  u  +  2v 


then 


9H 

at 


(VH) 


T  ZR 

at 


where 

3H 

3t 


R 


thus 


a2u 

ax2 


f  u  -  f  V 

y  z 


a  u  + 
u 


a  u  +  2  v 
21 


M— -  +  MBf  u  +  M$f  v 
3x2  y  z  . 


(3.21) 


By  the  curl  equations,  a  l2  =  a2i 


. 


> 
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Combining  this  with  (3.21)  and  (3.11)  yields 


V  = 


'  an“S+2MV!F  dX 

[  +  Mu0>  dx 


i  r 


u2  ( a  12  M3f  -  auf  )  +  v2(2M3f  -  a  12  f  ) 


dx 


u  v 


2 M 3 f  -  au  f  +a12(M3f_  -  f  ) 


dx 


(3.22) 


Here  again  use  can  be  made  of  the  inequality  employed  in 
chapter  2,  that  is 


±auv  £  |  a  |  |  u  v  |  £  |aj(u2  +  v2) 

2 

Using  (3.23)  and  integrating  by  parts,  (3.22)  becomes 


(3.23) 


i 


V  <  -  J  (M) 
J0  3x 

-  /'(— ) 

Jo  3x 


n 


al2  (M  +  1) 


2M  - 


a  (M  +  1) 
12  '  ' 


dx 

dx 


i 


/  u2Fu  dx  +  /  v2Fv  dx 

o  0 


(3.24) 


where 


Fu  =  a  12  M6fy  -  a „  fy  +  |2MBfy  -  au  fz  +  al2  (Mgfz  -  fy 

(3.25) 

Fv  =  2M6fz  -  a12fz  +  1 2  M  6  f  y  -  a„fz  +  au(HBfz  -  fy)| 

(3.26) 


Since 


u2  <  (1  -  x)  /( ill)*  dx 
“  Jo 


(3.27) 


(3.24)  can  be  wri tten  as 


> 
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V  < 


I 


(  — )  2  dx 
3x 


I 


(l-x)Fu  dx  -  a n 


a  ,,  ( M+ 1 ) 


2 


+ 


1 


1 


( 1 - x ) F  dx  -  2M 


■ai*AM+1^  (3.28) 

2  J 


Thus,  for  asymptotic  stability,  the  sufficient  conditions 
are 


(3.29) 


(3.30) 


Comparing  (3.29)  and  (3.30)  with  (2.35)  and  (2.36), 
respectively,  it  is  clear  that  both  sets  of  sufficient 
conditions  are  equivalent  if  the  following  relations  hold 

p  =  a  ,  p  =2,  p  =  a 

ri  n  2  ~  3  12 

Notice,  these  transformations  affect  Fu  and  Fv  as  well. 
There  is,  of  course,  no  loss  of  generality  in  having  set 
p  =2.  This  merely  acts  as  a  sort  of  scaling  factor,  in 


2 


that,  although  the  numbers  change  on  both  sides  of  the 
inequalities,  the  same  conclusions  hold  for  each  steady 
state.  It  amounts  to  multiplying  the  matrix  £  by  a 
constant.  The  results  are  shown  in  Table  4  which,  if  a 
comparison  is  made,  can  be  found  to  correspond  closely 
to  Table  1 . 


. 
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3.5  Another  Example:  The  Tubular  Reactor  Problem 

From  the  discussion  on  the  tubular  reactor 
problem  in  chapter  2,  recall  the  transformed  equations 
(2.54)  and  (2.55) 

ff  “  0  +  S(afy  -  £p)  +  anfz  (3.31) 

—  -  +  n(vf-  -  7P)  +  v5fv  (3.32) 

at  as2  z  4  y 

with  the  corresponding  boundary  conditions 

C ( t » 0 )  =  n(t,0)  =0, 

^(t.1)  -  -Jfp  e(t.l).  (3.33) 

as  2 

and  ~(t,l)  =  -^-p  n(t,l) 

as  2 

Proceeding  exactly  as  in  the  first  example  easily  leads  to 

an  equation  which  is  identical  to  (2.63)  except  that,  as 

in  the  previous  example,  the  following  changes  occur 

p  =  a  ,  p  =2  and  p  =  a 
1  n  2  3  12 

Thus,  it  is  obvious  that  the  same  Liapunov  functional 
will  result  if  the  same  inequalities  are  used.  The 
stability  results  obtained  using  the  above  parameters 
are  shown  in  Table  5. 


,1  ; 

\ 
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3.6  Conclusions 

The  method  has  been  applied  to  two  examples  in 

this  chapter.  In  both  cases  the  Liapunov  functional 

generated  was  the  same  (within  a  constant)  as  the  one  chosen 

for  use  in  the  generalized  quadratic  form  method.  This  was 

not  merely  by  coincidence,  but  rather,  as  planned.  It  was 

just  a  question  of  how  to  pick  the  values  of  a_.  Choosing 

other  values,  for  example,  in  either  case,  let  a  =0, 

12 

results  in  a  different  functional.  In  fact,  in  the  catalyst 
particle  problem  with  a12  =  0,  the  functional  which  will  be 
obtained  wi 1 1  correspond  to  a  constant  multiple  of  the 
functional  chosen  by  Berger  and  Lapidus.  Also  it  is  quite 
probable  that,  by  choosing  very  different  values  for  a  —  , 
more  complicated  types  of  Liapunov  functionals  will  result 
which  will  prove  stability  of  both  the  first  and  third 
steady  states  for  that  corresponding  system.  The  aim  of 
using  this  method  should  always  be  to  try  to  find  the 
simplest  form  of  Liapunov  functional  which  will  suffice. 

In  that  way,  the  job  of  maneuvering  the  equations  into  a 
form  from  which  sufficient  conditions  arise  is  very  much 
easier. 

When  using  the  variable  gradient  method  for 

i 

either  lumped  or  distributed  systems  it  is  wise  to  keep 
the  constant  and  variable  parts  of  the  a  ^  j  '  s  as  simple  as 
possible.  The  way  in  which  the  values  are  chosen  depends 


4 


' 

- 

» 

' 


very  much 
know! edge 


on  the  experience  of  the  analyst  and 
of  the  particular  plant  which  he  is 


h  i  s 

s  tudy i ng . 


I 
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CHAPTER  4 


CONCLUSIONS  AND  SUGGESTIONS 
FOR  FUTURE  RESEARCH 


4.1  Conclusions 

The  author  feels  that  there  are  two  main 
contributions  made  by  this  thesis.  Firstly,  the  use  of 
a  non-identity  matrix  P_  in  the  Liapunov  functional  and, 
secondly,  the  extension  of  the  variable  gradient  method 
to  distributed  parameter  systems. 

The  use  of  a  non-identity  matrix  in  the 
Liapunov  functional  is  a  direct  extension  of  the  quadratic 
form  used  frequently  as  a  Liapunov  function  for  determining 
stability  of  systems  described  by  ordinary  differential 
equations.  It  has  been  shown  here  that  there  are  advantages 
to  using  this  general  functional  form.  Whereas  the  special 
case  (P_=I_)  used  by  Berger  and  Lapidus  proved  stability  for 
steady  state  1  in  both  examples  treated,  no  conclusion 
about  stability  was  possible  for  the  other  two  steady 
states.  However,  by  using  a  non-identity  matrix  it  was 
possible  to  establish  stability  of  steady  state  3  as  well. 

It  has  also  been  shown  that  when  using  a  non-identity 
matrix  in  the  Liapunov  functional,  retention  of  all  the 
terms  in  the  inequalities  derived  is  necessary  as 
otherwise  inconclusive  results  are  obtained. 

The  extension  of  the  variable  gradient  method 
of  generating  Liapunov  functions  to  the  case  of 


' 


. 

■ 
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distributed  parameter  systems  has  been  described  through 
two  examples.  This  method  yielded  results  very  similar 
to  those  obtained  by  the  other  method  used  on  the  same 
examples.  It  has  been  shown  that  the  Liapunov  functionals 
which  were  arrived  at  by  the  variable  gradient  method 
are  very  similar  to  the  ones  used  in  the  first  method. 
This  is  due  to  the  choice  of  the  variable  parameters  at 
hand.  However,  the  end  results  were  obtained  somewhat 
more  systematically.  The  variable  gradient  method 
should  still  be  classified  as  a  trial  and  error  method 
although  each  trial  (if  more  than  one  is  necessary)  is 
carried  out  in  a  more  straightforward  manner. 

It  should  be  mentioned  that  in  both  of  the 
methods  presented  here,  there  is  the  requirement  that  the 
system  equations  be  at  least  of  second  order.  This  is 
necessary  for  using  a  matrix.  It  has,  however,  been  shown 
for  the  adiabatic  tubular  reactor  problem  that  results 
can  be  obtained  by  first  considering  the  more  general 
problem  and  then  reducing  the  final  inequalities  to  the 
adiabatic  case.  The  adiabatic  case  cannot  be  treated 
directly  by  either  method  since  the  original  transient 
equations  of  the  reactor  have  been  reduced  to  a  single 
equation. 

4 .2  Suggestions  for  Future  Research 

The  application  of  Liapunov  stability  theory 
to  distributed  parameter  systems  is  still  in  its  infancy. 


■ 


' 


' 


- 
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There  need  to  be  developed  general  methods  for  constructing 
Liapunov  functionals  for  various  classes  of  systems. 

Methods  which  could  deal  with  the  nonlinear  problem 
directly  without  linearization  would  be  more  reliable,  in 
that  they  would  probably  yield  sharper  results. 

The  variable  gradient  technique,  illustrated 
here  by  means  of  two  examples,  seems  promising  and 
should  be  further  developed  for  use  on  more  complicated 
types  of  systems. 
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APPENDIX  1 


APL  PROGRAMS 

For  the  interested  reader,  the  following  APL 

programs  are  listed  for  ready  reference.  It  was  mainly 

the  use  of  these  functions  that  yielded  the  results  of 

this  thesis. 

VRKGpjV 

V  R+T  RKG  X;TT  ;IT;J;Q1;Q2  ;Q3  ;S  N;XX;K  ;TSTEP 

[1]  ->(1  2  3  =p  ,RCON)/  3  4  5 

[2]  ->0  tR<-\0xpU+' INITIAL  TIME  MUS  T  B  E  SPECIFIED  IN  RCON} 

[3]  RCON+RCON , 1E~6 

[  4]  RCON+RCON 9lE~5x(tT+9T)Zp  9T"}-RCONll] 

[5]  S N+xRCO Nl  3]  xj T-I T+R CONI  l]x(r  =  T<-, T)lll 

[6]  J«-4  tR+x~TMI  N=TMIN+-0 . 02  xTl  p  T]  -IT 

[7]  Z-**(  4  ,  p  tX)  pX 

[8]  25TEP+SNxl/K  +  \(Tll]-IT)  ,0.02  xTlpTl-IT 

[9]  TT+IT 

[10]  Ql+TSTEPxTT  SET  XL  1;  ] 

[11]  Q2  +-TSTEP*  (  TT+TS  TEP*2  )  SET  Xll;l+Ql*2 

[12]  Q3<rTSTEPx(TT+ZSTEP*2  )  SET  Xll  il+Q2  *2 

[13]  X[/i>J[l;]-l*0.  1666666666666667x^1  +  (2  *Q2  +Q3)+TSTEP*( 
TT+TT+TS TEP)  SET  X[l;]+£3 

[14]  -*>(2  1  =eW- 1)/  16  17 

[15]  -*-9  t%5TEP*-TSTEP*2 

[16]  -*10,,X[1  3  ;  ]  <-Xl  3  1  ;] 

[17]  ->(  (XX+[  /  |  -/[  1]  Xl2  4  ;])<7?C0W[2  ])/22 +£[2  ]>£[l] 

[18]  XI  1  4  ;]«-*[  3  1  ;] 

[19]  TMIN+S  N*(  \TMIN)l  \TSTEP 

[2  0]  -►(  32  x  i  (  \TSTEP)  <\RCONL3l)  t2  1 ,  TS  TEP+TS  TEP*2 

[2  1]  +9  ,  J<-3 

[22  ]  TSTEP+SN*l/K  1  (  0  .  5xT[l~\-TT  )  >TS  TEP+TS  TEP+TS  TEP*XX< 
RCONl 2  ]*50 
[2  3]  TSTEP+-TSTEP*  2 
[2  4]  IT+TT 
[2  5]  4 

[2  6]  XZU1+X12  ;] 

[2  7]  -*(  ( | T[ 1 ] -T7) > |  TSTEPt2  )  /10 

[2  8]  R+R,  TTtXl2  ;  ] 

[2  9]  TSTEP^-SNxK  [  1]«-L  /X  [2  ]  ,  |  IT-T+l  iT 

[30]  +  ( 0*pT) /10-2*TSTEP=0 

[31]  ->0  ,  (  (pi?)*  l  +  (  pX)l2  ]  )  ,  l  +  (pX)[2  ]  )pR 

[32]  -*31,pD  +-'STEP  SIZE  HAS  B  ECO  bE  S  MALLER  THAN  INPUT  CON 
TROL  -  ’■RCONL  3  ]  • 


V 
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V  R+T  SET  X\N\B  \  A  \Q 

[ 1]  R+Xl 2 ] , ( 30x*C2 ])+22500000000000x(*Cl]- 

1. 4)x*-30vX[ 1] 

V 


VFEFNZUlV 


V  II+Z  FEFN  P\DFY  \DFZ 

Cl]  4«-2  .2  5x10*13 

[2  ]  DFY-*-*-3  0tZ[2  ;  ] 

[3]  DFZ+DFY  x3  0  x  ( l .  4-ZC2  ;])f0.  4xZ[2  ;]*2 

[4]  FE+(Pll]x(-A*DFY)-22  5  )  +  (P[3]x 

0.  *xA*DFY)+SS+0.  5x  |  (P[l]x^xPFZ)  +  (P[  3]x( 

0.  4x>l x  DFZ  )-  45  0 -Ax  DFY)+Pl2  ]x0.  H*A*DFY 

[5]  FN+(P  [2  ]x(  0.  H*A*DFZ)  -2  2  5  )  +  (  PC  3]  x  Ax  DFZ  )+SS 

C  6  ]  JT*-(  3  ,  p  Z  C  1;]  )  p  Z  C  1 ;  ]  , ( Z  C 1 ;  ]  *FE ) ,  Z  C 1 ;  ]  *FN 

V 


V INTEGRAL 1CD]V 


Cl]  "*3  x  \  2  -ppX 

C2  ]  -+0#pD«-,j?£QHn?ff  MATRIX  ARGUMENT' 
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VTEST [□]  V 


V 

Z  TEST  P;  Y 

Cl] 

II+Z  FEFN  P 

[2  ] 

Y+INTEGR  AL 1  II 

C3] 

-+(  CA+Yl  1]<P[  l]-2  x  |  P[  3  ]  )  / 5 

C4] 

’ CONDITION  A  FAILS' 

[5] 

+(CB+Y [2  ]  <P[2  ]-2  x  |  PC  3  ]  )/7 

C6] 

'CONDITION  B  FAILS' 

C7] 

+(  CA*CB  )  /10 

[8] 

'NO  CONCLUSION  LHS  (A)  =f;Y[l];» 

LHS  (B  )  =  »  ;Y[ 

2] 

[9] 

■>0 

[10] 

'THIS  SOLN.  IS  STABLE  ,  LHS  U)  =' 

r  ; Y[ 1 ]  ; '  LH 

S(B  )  =’  ;  Y[2  ] 
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APPENDIX  2 


DERIVATION  OF  EQUATIONS 
a)  Catalyst  Particle  Problem 

In  order  to  treat  a  one-dimensional  problem 
consider  the  catalyst  particle  to  have  the  form  of  a  flat 
slab  with  sealed  edges  as  shown. 


A  material  balance  over  this  element  yields 


|jRate  of  transport  i  nj  -  [Rate  of  transport  out] 

-  [Rate  of  reaction]  =  [Time  rate  of  change  ofj 

[concentration  I 


Ax  D 

3c 

3  C 

-  AxAAxi!l  =  iLLxAxAA 

( A2 . 1 ) 

3  A 

A+AA 

A 

d  x  3  x 

Dividing  through  by'  Ax  a  A  and  taking  the  limit  as  AA->0 
yields 

D  ii£.  -  ID.  =  i£  (A2.2) 

8  A  2  d  t  3  t 

where  D  =  the  effective  diffusion  coeffient 
c  =  concentration 

i  n 

—  =  per  unit  volume  reaction  rate 
d  x 


. 


. 

. 
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If  C  =  C0 

a  t  x  = 

L, 

( A2 . 

2) 

can  be  written 

32y 

J_ 

d  n  _ 

l 

ay 

3  X  2 

Dc0 

dt 

D 

3  T 

where  y 

ii 

n 

o 

o 

This  may 

be  normalized  by 

letting  x  =  x  /  L . 

Then 

3  2y 

ki 

d  n  _ 

.  £X 

3  x  2 

Dc0 

d  t 

D 

3  T 

From  the  work  of  Weisz  and  Hicks  i 


that 


=  k  c  y  exp 
dx  oo" 


y(z-1) 

z 


where 


k  =  rate  constant  at  the  boundary 
o 


z  =  1/1 


o 


T  =  temperature 


T  =  Tq  at  the  boundary 


combining  (A2.4)  and  (A2.5)  yields 


a2y  _  L2k0 


o  r 


where 


and 


3x2 

3  2y 

3X2 


D 


y  exp 


y 


(z-1) 


=  3y 

D  3  t 


2  - 


-  <pzy  exp 

i  ^  k 

L  Ko  


Y  (  z-1  ) 


D 


=  iX 

3  t 


Thiele  modulus 


t  =  kx 
L2 


Equation  (A2.7)  is  the  same  as  equation  (2.3) 


(A2.3) 

(A2.4) 

t  follows 

(A2.5) 


(A2.6) 

(A2.7) 


V 


t 
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A  similar  temperature  balance  on  the  element 

yields 


A*K 

3T 

-  il 

-  AxAXxHxdn  = 

ilpc  xAxAX 

(A2.8) 

3  X 

X  +  AX 

3  X 

X 

d  t 

3  T  P 

where  K  =  effective  thermal  conductivity 

H  =  molar  heat  of  reaction 
p  =  density 
0^  =  heat  capacity 


In  the  limit  (A2.8)  can  be  written 


K—  -  =  pc  —  (A2.9) 

3  X 2  d  t  3  t 


normalizing  yields 


32z 

L2H 

dn  _  L2pc 

—  -  - ~p 

3  Z 
— 

(A2.10) 

3x2 

KT 

d  t  K 

3  T 

0 

or 

Z'zz 

L2  H 

k  c  v  exp 

0  0^ 

Y ( Z-  1  ) 

_  L2P  C  3  z 

(A2. 11) 

3X2 

KT„ 

Z 

K  P  3t 

which  can  be  simplified  to 


where 


32z 

3x2 


+  <j>23  y 


exp 


y ( z-  1 ) 

z 


3 

N 


c0HD 

KTo 

PcpD 

K 


-ft1) 

o  max 

=  Lewis  number 


(A2. 12) 


Equation  (A2.12)  is  the  same  as  (2.4) 
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b)  Adiabatic  Tubular  Reactor 

A  similar  derivation  for  equations  (2.40)  and 
(2.41)  can  be  found  i  n  ^ .  For  clarity  some  of  the  constants 
are  defined  here. 


y  =  c/c0 
z  =  T/T0 


where  c0,T0  occur  at  the  inlet  of  the  reactor. 


O’  0 

• 

N 

=  vl 

P 

D 

where 

v  = 

fluid  velocity  within  reactor 

1  = 

reactor  length 

a  = 

1  2  K 

D 

Q  = 

E/RTo 

where 

E  = 

activation  energy 

and 

R  = 

gas  constant 

v  = 


B  a 

_  -He 


o 


TqP  c. 


where 


B 


■ 


■ 


